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Abstract 



We outline basic principles of canonical formalism for the Nambu mechanics — a 
generalization of Hamiltonian mechanics proposed by Yoichiro Nambu in 1973. It 



is based on the notion of Nambu bracket, which generalizes the Poisson bracket— 
a "binary" operation on classical observables on the phase space, to the "mul- 
tiple" operation of higher order n > 3. Nambu dynamics is described by the 
phase flow given by Nambu-Hamilton equations of motion — a system of ODE's 
which involves n — 1 "Hamiltonians" . We introduce the fundamental identity 
for the Nambu bracket — a generalization of the Jacobi identity, as a consistency 
condition for the dynamics. We show that Nambu bracket structure defines an 
hierarchy of infinite families of "subordinated" structures of lower order, includ- 
ing Poisson bracket structure, which satisfy certain matching conditions. The 
notion of Nambu bracket enables to define Nambu-Poisson manifolds — phase 
spaces for the Nambu mechanics, which turn out to be more 'rigid" than Pois- 
son manifolds — phase spaces for the Hamiltonian mechanics. We introduce the 
analog of the action form and the action principle for the Nambu mechanics. In 
its formulation dynamics of loops (n — 2-dimensional chains for the general n-ary 
case) naturally appears. We discuss several approaches to the quantization of 
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Nambu mechanics, based on the deformation theory, path integral formulation 
and on Nambu-Heisenberg "commutation" relations. In the latter formalism we 
present explicit representation of Nambu-Heisenberg relation in the n = 3 case. 
We emphasize the role ternary and higher order algebraic operations, and math- 
ematical structures related to them, play in passing from Hamilton's to Nambu's 
dynamical picture. 

1 Introduction 

In 1973 Nambu proposed a profound generalization of classical Hamiltonian mechanics [|TJ. 
In his formulation a triple (or, more generally, n-tuple) of "canonical" variables replaces a 
canonically conjugated pair in the Hamiltonian formalism and ternary (or, more generally, 
n-ary) operation — the Nambu bracket — replaces the usual Poisson bracket. Dynamics, ac- 
cording to Nambu, is determined by Nambu-Hamilton equations of motion, which use two 
(or, more generally, n — 1) "Hamiltonians" and replace canonical Hamilton equations. Corre- 
sponding phase flow preserves the phase volume so that the analog of the Liouville theorem 
is still valid, which is fundamental for the formulation of statistical mechanics [|l|. 

Nambu's proposal was partially analyzed in papers || [|. In [0] it was shown that 
particular example of Nambu mechanics can be treated as a six-dimensional degenerate 
Hamiltonian system with three constraints and a Lagrangian being linear in velocities. In 
|| it was shown that one can use a four-dimensional phase space as well. However, until 
recently, there were no attempts to formulate the basic principles of Nambu mechanics in 
the invariant geometrical form similar to that of Hamiltonian mechanics ||. In this paper 
we develop the basics of such formalism and display novel mathematical structures which 
might have physical significance. 

We start by formulating the fundamental identity (FI) for the Nambu bracket as a con- 
sistency condition for the Nambu's dynamics. As a corollary, it yields the analog of the 
Poisson theorem that Poisson bracket of integrals of motion is again an integral of motion. 
Based on FI, we introduce the notion of Nambu-Poisson manifolds, which play the same role 
in Nambu mechanics that Poisson manifolds play in Hamiltonian mechanics. 

We show that Nambu bracket structure contains an infinite family of "subordinated" 
Nambu structures of lower degree, including Poisson bracket structure, with certain matching 
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conditions. This implies that Nambu bracket structure is, in a certain sense, more "rigid" 
then the Poisson bracket structure. It can be seen explicitly by showing that FI imposes 
rather strong constraints on possible forms of Nambu bracket. In addition to quadratic 
differential equations, it also introduces an overdetermined system of quadratic algebraic 
equations for the Nambu bracket tensor. This is the novel feature in comparison with 
the Poisson bracket case, where we have differential constraints only. Additional algebraic 
requirements are "responsible" for the above-mentioned "rigidity" of the Nambu bracket 
structure in comparison with the Poisson bracket structure. Specifically, whereas any skew- 
symmetric constant 2-tensor yields a Poisson bracket, this is no longer true for the Nambu 
bracket. This manifests rather special nature of Nambu mechanics. From our point of view 
such exclusiveness should be considered as an advantage of the theory rather then of the 
contrary. 

Next, we develop canonical formalism for the Nambu mechanics. It is based on the analog 
of Poincare-Cartan integral invariant — action form, which is a differential form of degree 
n — 1 on the extended phase space. It has properties similar to those of the usual Poincare- 
Cartan invariant and enables us to formulate principle of least action for Nambu mechanics. 
Instead of all possible "histories" of the Hamiltonian system considered in action principle 
for classical mechanics (all paths connecting initial and final points in the configuration 
space), in Nambu mechanics one should consider all n — 1-chains in the extended phase 
space whose "time-slices" are closed n — 2-chains satisfying certain boundary conditions. 
We define classical action as an integral of generalized Poincare-Cartan invariant over such 
n — 1-chains and prove the principle of least action: "world- sheets" of a given n — 2-chains 
under the Nambu-Hamilton phase flow are extremals of the action. 

Finally, we briefly discuss quantization of Nambu mechanics. This problem, first consid- 
ered by Nambu [Q, is still outstanding. We indicate several possible approaches towards its 
solution. In particular, we construct a special representation of Nambu-Heisenberg "com- 
mutation" relations, which is similar in spirit to the representation of canonical Heisenberg 
commutation relations by creation-annihilation operators in the space of states of the har- 
monic oscillator. In our realization states are parametrized by a lattice of algebraic integers 
in a cyclotomic field for the cubic root of unity, whereas states of the harmonic oscillator are 
parametrized by non-negative rational integers. 

Now let us explain the main ideas of the paper in more detail. We start with the 
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simplest phase space for Hamiltonian mechanics — a two-dimensional Euclidean space 1R 2 
with coordinates x, y and canonical Poisson bracket 

' dx dy dy dx d{x, y) 

This bracket satisfies Jacobi identity 

{/i, {/ 2 , /s}} + {/a, {fx, / 2 }} + {/a, {/s, A}} = 0, 
and gives rise to a Hamilton's equations of motion 

s = <*■'>• 

where / is a classical observable — a smooth function on the phase space — and if is a Hamil- 
tonian. 

Generalization of this example leads to a concept of Poisson manifolds — smooth mani- 
folds endowned with a Poisson bracket structure satisfying skew-symmetry condition, Leibniz 
rule and Jacobi identity. 

Canonical Nambu bracket |IJ is defined for a triple of classical observables on the three- 
dimensional phase space 1R 3 with coordinates x, y, z by the following beautiful formula 

{ f f f 1 _ 9(fi, /2, /3J 

\Ji,j2, J3} — —^7 \— , 

d{x,y,z) 

where the right-hand side stands for the Jacobian of the mapping / = f 2 , fs) : R 3 1— > ]R 3 . 
This formula naturally generalizes usual Poisson bracket from binary to ternary operation 
on classical observables.Q Generalized Nambu-Hamilton equations of motion involve two 
"Hamiltonians" Hi and H2 and have the form 

Corresponding phase flow on the phase space is divergence-free and preserves the standard 
volume form dx A dy A dz - analog of the Liouville theorem for Nambu mechanics |l] . 

In Section 2 we prove that canonical Nambu bracket satisfies the following fundamental 
identity (and its generalizations for the n-ary case)0 

{{fli fl, fz\i fii h} + {/3> {fli /2> ^4}, fs} + {/3> fx, {/l, /2> A}} = /2; {^3; ^4, /s}}- 



1 M. Flato §| informed me that, apparently, Nambu introduced this bracket in order to develop a "toy 
model" for quarks considered as triples. 

2 This relation was also independently introduced by F. Bayen and M. Flato M. 



This formula might have been considered as the most natural (at least from a "dynamical" 
point of view) generalization of the Jacobi identity. It yields (see Theorem 3) an analog of 
classical Poisson theorem that Poisson bracket of two integrals of motion is again an integral 
of motion. Generalized version of FI for the n-ry case (where Jacobian of the mapping / = 
(fi, . . . , f n ) : ]R n i— > R n defines canonical Nambu bracket of order n) enables us to introduce 
Nambu-Poisson manifolds as a smooth manifolds with Nambu bracket structure — an n-ary 
operation on classical observables satisfying skew-symmetry condition, Leibniz rule and FI. 
We show that by fixing some of the arguments in the Nambu bracket of order n one can 
get brackets of lower order which still satisfy fundamental identity (i.e. are Nambu brackets) 
and, in addition, satisfy certain matching conditions for different choices of fixed arguments. 
This hierarchical structure of Nambu bracket shows that it is more "rigid" concept yen that 
of Poisson bracket. We discuss explicit conditions FI imposes on corresponding n-tensor 
of the Nambu bracket. We show, contrary to the Poisson case, that algebraic part of FI 
substantially reduces possible Nambu structures with constant skew-symmetric n-tensor. 
Specifically, we show that bracket on the phase space B? n = ©™ =1 R 3 , defined in [§ (and 
used in H) as a direct sum of canonical Nambu brackets on IR 3 , does not satisfy FI and, 
therefore, is not a Nambu bracket. Next, we discuss linear Nambu brackets. We show that 
they naturally lead to a new notion of Nambu-Lie "gebras" , which generalizes Lie algebras 
for the n-ary case. We close Section 2 by presenting several simple examples of evolution 
equations which admit Nambu formulation. 

In Section 3 we introduce the analog of Poincare-Cartan integral invariant for Nambu 
mechanics — a differential form of degree n — 1 on the extended phase space. In the simplest 
three-dimensional example described above it is given by the following 2-form on IR 4 

w ( 2 ) = xdy Adz - H 1 dH 2 A dt. 

We prove in Theorem 6 that the vector field of Nambu-Hamilton phase flow in the extended 
phase space is the line field of the 3- form duj^ 2 \ so that integral curves are its characteristics. 
We define a classical action as an integral of Poincare-Cartan action form over n — 1-chains, 
and in Theorem 7 prove the principle of least action. It states that n — 1-chains — "tubes" 
of integral curves of Nambu-Hamilton phase flow "passing through" a given n — 2-chains — 
are the extremals of the action. These results generalize the basic facts lying the foundation 
of Hamiltonian mechanics (see, e.g., ||). 
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In Section 4 we discuss possible approaches to the quantization of Nambu mechanics. 
Though we mention those based on the deformation theory and Feynmann path integral, 
our main result is the explicit construction of a special representation of Nambu- Heisenberg 
commutation relation. For the ternary case this relation has the form (see JTJ) 

[A l} A 2 , A 3 ] = A X A 2 A 3 - A X A 3 A 2 + A 3 A X A 2 - A 3 A 2 A 1 + A 2 A 3 A l - A 2 A X A 3 = cl, 

where A\,A 2 ,A 3 are linear operators, J is a unit operator and c is a constant. Let Z[p] be 
a lattice of algebraic integers in quadratic number field Q[p], where p 3 = 1. We prove in 
Theorem 8 that Nambu-Heisenberg relation can be represented by operators A± acting in 
the linear space H = {\co), ou G Z[p]} by the following simple formula 

Ai\u) = (w + 1 + p)\u + 1), A 2 \u) = (u + p)\u + p), A 3 \u) = lj\uj + p 2 ). 

This realization should be compared with canonical representation of Heisenberg commu- 
tation relations in the space of states of harmonic oscillator given by creation-annihilation 
operators. Contrary to the latter case, our representation does not have a vacuum vector 
(at least in the conventional sense) and we do not know whether the analog of Stone-von 
Neumann theorem (unitary equivalence of irreducible representations) is still true. 

Acknowledgments. I appreciate stimulating conversations and discussions with Moshe 
Flato, Igor Krichever, Nicolai Reshetikhin, Han Sah and Daniel Sternheimer. I am grateful to 
Larry Lambe for performing symbolic computations in connection with a system of quadratic 
algebraic equations described in Section 2. I should also mention that Theorem 8 was 
obtained as early as in 1977. Its publication was deliberately postponed since at that time 
it did not "fit" into a general picture. 

This work was partially supported by the Paul and Gabriella Rosenbaum Foundation. 

2 Nambu-Poisson Manifolds 

In Hamiltonian mechanics a smooth manifold X is called a Poisson manifold and its function 
ring A = C°°(X) — algebra of observables, if there exists a map { , } : A <g> A 1— > A with the 
following properties. 
1. Skew-symmetry 

{/i,/a} = -{/*, A}, 
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for all /i, h e A. 

2. Leibniz rule (derivation property) 



{/l/2,/s} = /l{/2,/s} + /2{/l,/s}- 



3. Jacobi identity 



{fi, {/ 2 , / 3 }} + {/ 3 , {A, / 2 }} + {/ 2 , {/ 3 , A}} = o, 



for all /i,/ 2 ,/ 3 e A 



Corresponding "binary operation" { , } on A is called Poisson bracket and plays a 
fundamental role in classical mechanics. Namely, according to Hamilton, dynamics on the 
phase space X is determined by a distinguished function H e A — a Hamiltonian, and is 
described by Hamilton's equations of motion 



When solution to Hamilton's equations exists for all times t E R and all initial data (this 
is so when X is compact), it defines Hamilton phase flow x ^ g t {x), x G X, and evolution 
operator U t : A i— > A, 



Hamilton's dynamical picture is consistent if and only if evolution operator U t is an 
isomorphism of algebra of observables A. This means that U t is an algebra isomorphism, 
i.e. Utififz) = Ut{fi)Ut{fo) and, in addition, preserves the Poisson structure on A, i.e. 
Ut({fi, / 2 }) = {Ut(fi), ^t(/ 2 )}- It is easy to see (using the standard uniqueness theorem 
for ODE's) that the first property is equivalent to the Leibniz rule and the second one is 
equivalent to the Jacobi identity. We summarize these well known results as the following 
theorem. 

Theorem 1 Evolution operator in Hamilton's dynamical picture is an isomorphism of 
algebra of observables A = C°°(X) if and only if the phase space X is a Poisson manifold. 

The basic examples of Poisson manifolds are given by two-dimensional phase space X = 
R 2 with coordinates x, y and Poisson bracket 



df 



{HJ}, feA. 



dt 



U t (f)(x) = f(g t (x)), xeX, feA. 



{/l,/ 2 } 
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and by its generalization — X = R 27V with coordinates x\,... , XN,yi, ■ ■ ■ ,Vn and Poisson 
bracket 

U 2 ~i dxi dyi dyi dxi 
Geometrically, Poisson manifold X is characterized by a Poisson tensor r\ - a section of the 
exterior square A 2 TX of a tangent bundle TX of X, which defines Poisson structure by the 
formula 

{fi,f2} = v(dfi,df 2 ). 

Jacobi identity is equivalent to the property that n has a vanishing Schouten bracket with 
itself. In local coordinates (xi, . . . ,xn) on X Poisson tensor n is given by 

N d d 

and Jacobi identity identity takes the form 

#v drj jk dr/ki %^ 

S ( ""a^ + ^&r + ""&7 ) = ' 

for all i,j,k = \,...,N. 

Dynamics according to Nambu consists in replacing Poisson bracket by a ternary (n- 
ary) operation on algebra of observables A and requires two (n — 1) "Hamiltonians" Hi, H 2 
(Hi, . . . , H n -i) to describe the evolution. This dynamical picture is consistent if and only if 
evolution operator is an isomorphism of algebra of observables. Therefore, we propose the 
following definition. 

Definition 1 Manifold X is called a Nambu- Poisson manifold of order n if there exists 
a map {,...,} : A®" i— > A — generalized Nambu bracket of order n, satisfying the following 
properties. 

1. Skew-symmetry 

{fl, ■ ■ - ,fn} = (-l) e(<j) {/a(l), • • • , fa(n)}, 

for all fi, . . . , f n G A and a G Symm{n), where Symm(n) is a symmetric group of n elements 
and e(<r) is the parity of a permutation a. 

2. Leibniz rule 



{/l/2, /3, • • • , /n+l} — /l{/2, /3, • • • , /n+l} + .M/l, ^3, • • • , /n+l}, 

8 



for all /i, ...,/„ + i G A. 

3. Fundamental identity (FI) 

{{/l, • • • , /n-l, /n}, /n+1, • • • , /2n-l} + {/n, {/l, • • • , /n-l, /n+l}, /n+2, • • • 5 /2n-l} 
+ • • • + {/n, • • • 5 /2n-2, {/l, • • • , /n-l) /2n-l}} = {/l, • • • , /n-l, {/n, • • • , /2n-l}}, 

for all /1, ...,f 2 n-l e A 

Remark 1 Nambu bracket structure of order n on a phase space X induces infinite 
family of "subordinated" Nambu structures of orders n — 1 and lower, including the family 
of Poisson structures. Indeed, consider n = 3 case and for if G A define the bracket { , }h 
on X as 

{^0}^ = {ii,^,0} 

for all ip,(/> E A. Setting in FI fi = fs = H, we see that it turns into the Jacobi identity 
for the family of brackets { , }h parametrized by observables H. Conversely, such family 
of Poisson brackets gives rise to Nambu bracket if certain matching conditions are satisfied. 
Namely, for any e A define 

D?(f) = {4>,f}H, feA, 

which is the derivation of the Poisson bracket {, , }h (Jacobi identity). Then the family 
{ , } H , H G A of Poisson brackets on X gives rize to a Nambu bracket, defined as 

{/l, /2, fz} = {/2, /3}/i, 

if and only if 

for all Hi,H 2 G A. Indeed, it is easy to see that this equation — a derivation property of 
D^ 1 with respect to the whole family of Poisson brackets { , }h — is equivalent to FI (if one 
identifies f\ = Hi, f 2 = <p, fz = ip, fi = x an d /s = H 2 ). Moreover, this condition for the 
case Hi — H 2 — H is equivalent to the Jacobi identity for the bracket { , The same is 
true for the general case, where for all Hi, . . . H n _ k G A the assignment 

{/i, • • • , fk}n 1 ...H n _ k = {Hi, . . . , H n ^ k , fi, . . . fk) 
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defines a hierarchy of subordinated Nambu structures of orders k = 2, . . . , n—1, parametrized 
by the elements in A n ~ k A. They all satisfy FI (which follows from FI for the "basic" structure 
of order n) and matching conditions of the same type as above. 

Dynamics on Nambu-Poisson manifold is determined by n — 1 functions Hi, . . . ,H n _i 
and is described by generalized Nambu-Hamilton equations of motion 

f t ={H h ...,H n „i,f}, feA. (1) 

Corresponding Nambu-Hamilton phase flow g t defines evolution operator Ut, Ut(f)(x) = 
ftfx), x e X, for all feA. 

The following theorem clarifies "dynamical" meaning of a concept of Nambu-Poisson 
manifolds. 

Theorem 2 Evolution operator in Nambu 's dynamical picture is an isomorphism of 
algebra of observables A = C°°(X) if and only if the phase space X is a Nambu-Poisson 
manifold. 

Proof. We need to prove that 

U t ({h,...J n }) = {U t (f 1 ),...,U t (f n )}. (2) 

Since (Q) is obviously valid at time t = 0, it sufficient to show that it both sides satisfy 
the same evolution differential equation. Denoting by L e Vect(X) the vector field of the 
Nambu-Hamilton flow <?*, i.e. 

L(f) = {Hi,...,H n _i,f}, 
we can express i-derivative of (0) as 

L({fl, ...,fn}) = {L(h), f 2 , . . . , f n } + {fi, L(f 2 ), ...,f n } + {fl, f 2 ,..., L(f n )}, 

which is nothing but FI specialized for the functions Hi, ... , H n -i, fi, ■ ■ ■ , f n - 1=1 

Remark 2 Equivalent formulation of FI used in the previous proof can be stated ex- 
plicitly that for any elements Hi, ... , H n _i the mapping L : A i— > A is a derivation of the 
Nambu bracket. 

Definition 2 Observable F e A is called integral of motion for Nambu-Hamilton system 
with Hamiltonians Hi, . . . , H n _i, if 

{Hi,...,H n „i,F} = 0. 
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As an obvious corollary of FI we get the following result. 

Theorem 3 Nambu bracket of n integrals of motion is an integral of motion. 

Remark 3 As we have seen in Remark 1, Nambu bracket structure of order n contains an 
infinite family of subordinated structures of lower degree, including Poisson bracket structure, 
with matching conditions between them. Therefore one might suspect that Nambu structure 
should be more "rigid" then its Poisson counterpart. This can be seen by comparing FI of 
order n with Jacobi identity — its special case when n — 2. As we know, the left hand 
side of Jacobi identity, considered as a map from A ® A ® A into A, is a deriviation with 
respect to every argument. However, for the general case n > 3 the difference between 
the left and right hand sides of FI, considered as a map from 2n — 1-fold tensor product 
A <g> . . . (gi A into A, is a derivation only with respect to the arguments f2 n -i and not 

to f\, . . . , fn-i- This is because these two groups appear in FI in a different way. Namely, 
analyzing FI, it is easy to see that observables from the first group appear only twice under 
the "double Nambu bracket", whereas there n such terms for every member in the second 
group^. Consider, for instance, such terms for the argument /^n-i- Denoting by L ix ^ in _ x 
vector field corresponding to Nambu-Hamilton phase flow with Hamiltonians . . . , fi n _ l , 
we can arrange these terms as a commutator \L n ...in-i-> -^i...n-i](/2n-i), so they are still given 
by a vector field action. This proves the derivation property with respect to the argument 
fm-\ (cf. standard arguments in Hamiltonian mechanics Same arguments apply to all 
members of the second group, but not to those of the first group. As we shall see below, this 
feature of FI implies strong constraints on the possible forms of Nambu bracket. 

Geometrically, Nambu structure of order n can be realized as 

{fu •••,/«} = v(dfi,---,df n ), 

where i] is a section of the n-fold exterior power A n TX of a tangent bundle TX. In local 
coordinates (x\, . . . , xn) on X Nambu tensor rj is given by 

N d d 

V= E Vh...i n (x)-K— A... A — 

h,...,in=l 11 tn 

and should satisfy FI. As we mentioned earlier, FI implies strong constraints on n-tensor rj. 

3 This is not surprising since, according to interpretation in Remark 2, members from the second group 
may be considered as a Hamiltonians for Nambu-Hamilton equations of motion, whereas members from the 
first group are "just" observables. 



11 



First, taking into the account Remark 3, we see that all terms containing second deriva- 
tives of fi, . . . , / n _i should vanish. This results in the following system of quadratic algebraic 
equations 

% + P(A% = 0, (3) 
for all multi-indices i = . . . , i n } and j = {ji, . . . ,j n } from the set {1, . . . , N}, where 

Nij — Vil—i n Vjl—jn r lj n i\iz---in r lji-..jn-li2 + • • • + Vjni2~-in-lilVjl~-jn-lin ~~ Vjni2--inVjl—jn-lh- (^) 

and P is a permutation operator which interchanges first and n + 1-th indices (i.e. i\ and 
ji) of a 2n-tensor N. 

Second, all terms containing first derivatives of fx, . . . , fin-\ must vanish. This yields the 
following system of quadratic differential equations 

El ®Vjl...jn , 'I'll -j- :'.' , , dVjl-jn-lin \ _ r. /r\ 

\ilh 2 ...i„ T 'lj n h S ...l n T ■ ■ ■ T '/jn«2...«n-l« Q J U 5 W 

for all indices Z2, . . . , 2 n; Ji ; ■ ■ ■ ijn = 1> • • • j 

iV. 

Thus skew-symmetric n-tensor 77 defines Nambu bracket of order n if and only if it satisfies 
equations (H) - (!)■ 

This shows significant difference between Nambu and Hamiltonian formulations — a con- 
stant skew-symmetric n-tensor rj (which obviously satisfies (|5])) for n > 3 no longer "auto- 
matically" defines Nambu bracket! To do so, it must satisfy algebraic constraints ([]) - (Q). 
One may wonder whether there are any solutions at all to the algebraic-differential system 

(I) - (D? 

The following geometric interpretation of tensor N provides the simplest examples. Let V 
be the iV-dimensional linear space and V* be its dual space. Any constant skew-symmetric 
n-tensor rj can be interpreted as an element in the linear space A n V, which we also will 
denote by 77. With every 77 G A n V one can associate a map N : A n_1 V* 1— > A n+1 V, defined 
by the formula 

Na = ijrf) A 77 G A n+ V, 

for all a G A n ~ 1 V*. Here i a (rf) G V is given by (i a (r]),v*) = (77, aAt> *) for any ?j* G and ( , ) 
stands for the pairing between V and V*. It is well known that equation N = is equivalent 
to the condition that element 77 G A™^ is decomposable, i.e. there exist v 1, . . . v n G V such 



12 



that 7] = v\ A . . . A v n , and it is easy to verify that in coordinates the map N is represented 
by the 2n-tensor N given by formula (f|) . 
Thus we have proved the following result. 

Theorem 4 Let V be a linear space. Any decomposable element in A n V endows V with 
the structure of a Nambu-Poisson manifold of order n. 

In particular, consider Nambu's original example , when X = R n is a phase space with 
coordinates x\, . . . , x n and "canonical" Nambu bracket is given by 

( t t l d(fi, . . . , f n ) , . 

{/l "--' /n} = ^ — x~y (6) 

where the right-hand-side stands for the Jacobian of the mapping / = . . . , f n ) : R n \— > 
R n . From Theorem 4 we get 

Corollary Euclidean space R n with canonical Nambu bracket of order n is a Nambu- 
Poisson manifold. 

Proof. Canonical Nambu bracket of order n is given by totally anti-symmetric n-tensor 
Vh-in = e ii...i„i which corresponds to the volume element of R n and, therefore, is decompos- 
able. □ 

Remark 4 We define canonical transformations as a diffeomorphisms of the phase space 
which preserve Nambu bracket structure. For the Nambu's example linear canonical trans- 
formations form a group SL(n, R). 

Remark 5 In [[L|] Nambu also considered another example — "direct sum" of canonical 
brackets. In this case X = R 3n = ®^ =1 R 3 and the bracket is given by the following element 



7] = ei A e 2 A e 3 + e 4 A e 5 A e 6 + . . . + e 3n _ 2 A e 3n _i A e 3n G A R' 



: >F> : >" 



where ej, i — 1, . . . , 3n is a basis in R 3 ™ induced by a standard basis in R 3 . (This bracket 
was used in Q to write equations of motion of a particle interacting with SU(2) monopole). 
However, it is easy to see that such tensor r] does not satisfy system (|3]) and, therefore, 
does not define Nambu bracket! This "explains" Nambu's observation JTJ that linear canon- 
ical transformations for this bracket "decouple", i.e. form a direct product of n copies of 
SL(3, R) — a fact he considered to be rather disappointing. 

To summarize this discussion, we see that there is significant difference between Nambu 
and Hamiltonian mechanics with Nambu formulation being "more rigid" . 
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The problem of constructing other examples of Nambu-Poisson manifolds is of great 
importance. Not trying to address it here (but hoping to to return to it later on), we will 
mention a just a few. 

Remark 6 It seems that in the constant case there should be other examples of Nambu 
brackets besides those given by decomposable tensors, since equations @ do not immediately 
imply that iV = 00. 

Another class of examples is provided by a non-constant tensors rj. As we know, linear 
Poisson bracket structure (Poisson structure on linear space such that Poisson bracket of 
linear functions is again linear) is equivalent to a Lie algebra structure on the dual space. 
One can ask what kind of structure does linear Nambu bracket introduce? 

Definition 3 A vector space V is called Nambu-Lie "gebra" of order n if there exists a 
map — Nambu-Lie bracket — [ .,...,.] : A n V i— > V such that 

[[fl, . . .,V n -i,V n ],V n+ x, . . .,V 2n -l] + K, [Vl, ■ • ■,V n -i,V n+ x],V n+ 2, ■ ■ .,V 2n -l] 
+ • • • + K, . . . ,t>2n-2, [Vx, ■ ■ ■ ,V n -X,V 2n -x]] = fal, • • • ,V n -X, K, • • • ,V 2n -l]], 

for all vi, . . . ,v 2n -i e V. 

Remark 7 Definition of Nambu-Lie "gebras" can be stated in equivalent form as a 
condition that for any Vi, . . . , t> n _i G V "adjoint" map [vx, ■ ■ ■ , f ra -i, ] : V i— > V is a derivation 
with respect to the Nambu-Lie bracket [,..., ]. 

Theorem 5 Linear Nambu structures of order n are in one-to-one correspondence with 
Nambu-Lie "gebras" of order n on the dual space. 

Proof. Fundamental identity for linear functions is nothing but FI in the definition of 
Nambu-Lie "gebras". □ 

In coordinates, linear Nambu structure of order n is given by n + 1-tensor c^ „■ : 

" 41 •••on 

N 

{ x ili • • • ? x i n } ~ C i 1 ...i n X k- 
k=X 

"Structure constants" c\ ' in satisfy overdetermined system of quadratic algebraic equations 
which follow from (^|) — (|5|). We are planning to analyze this interesting structure elsewhere. 

Larry Lambe, using symbolic computations technique, present an explicit form of a Groebner basis of a 
polynomial ideal of all relations (||) in the case n — 3, N = 6. It follows from his analysis that in this case 
all solutions are represented by decomposable tensors. 
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Now we present several simple examples of dynamical systems which admit Nambu for- 
mulation. 

Example 1 It goes back to Nambu |l[ and is given by Euler equations for the angular 
momentum of a rigid body in three dimensions. These equations admit both Hamiltonian 
formulation with respect to the linear Poisson bracket on R 3 = su(2)*, where Hamiltonian 
is given by kinetic energy and Nambu formulation with respect to the canonical ternary 
Nambu bracket on 1R 3 and two Hamiltonians - kinetic energy and total angular momentum. 

Example 2 Lagrange system (sometimes called Nahm's system in the theory of static 
5't/(2)-monopoles, see, e.g., [§, |10| ) on R 3 , which is given by the following equations of 
motion 

dx\ dx 2 dxs 

-7T = X 2 X 3 , — = X1X3, — = X±X2, 

at at at 

can be written in Nambu form 

dx ■ 

= {H u H 2 ,x i } ) i = 1,2,3, 

where H\ — x\ — x\, H 2 = x\ — x\. Integrals H\ and H 2 confine the phase flow to the 
intersection of two quadrics in R 3 - a locus of an elliptic curve so that the system can be 
integrated by elliptic functions. There exists another system with quadratic nonlinearity — 
so-called Halphen system, which is related to the Lagrange system (see, e.g., |9|, |I0|). It can 
be integrated in terms of modular forms (see, e.g., |T(|) and does not admit global (single- 
valued) integrals of motion. However, this system has two multi-valued integrals which play 



the role of two Hamiltonians in Nambu formulation ||11|| . 

Example 3 Let X = R n be the phase space with canonical Nambu bracket of order n 
and let the elementary symmetric functions of n variables Xi, . . . , x n to be the Hamiltonians: 



Hi — si, . . . , H n _i — s n _ 



1: 



where (x — x\) . . .{x — x n ) = x n — s\x n . . . ± s n . Nambu-Hamilton equations of motion 

^ = {H 1 ,...,H n - 1 ,x i }, i = l,...,n, 



can be written as 



dx, d n f . 
— — = — — , 1 = 1, ...,n, 
dt oxf 
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where 



l<i<j<n 



which might suggest a generalization of a gradient flows. 

3 Canonical Formalism and Action Principle 

Here we extend canonical formalism of Hamiltonian mechanics, based on Poincare-Cartan 
integral invariant and on the principle of least action (see, e.g., ||), to the case of Nambu 
mechanics. We illustrate essential features of our approach on the simplest example of 
Nambu- Poisson manifold X = 1R 3 with canonical Nambu bracket. 

Let X = 1R 4 be the extended phase space with coordinates x, y, z, t. 

Definition 4 The following 2- form to^ on X 



uj {2) = xdy Adz - H 1 dH 2 A dt, 



(7) 



is called generalized Poincare-Cartan integral invariant — action form for Nambu mechanics 
(cf. §). 

Consider Nambu-Hamilton equations with Hamiltonians Hi and H 2 

df 



where 



and 



Denote by 



f t = {H u H 2 J} = L(f), 



d d d 

L = Li — + L 2 — + L 3 — e Vect(X), 
ox ay oz 



T d(Hi,H 2 ) _d(H u H 2 ) _ d(Hi,H 2 ) 

L l — — K7 \ — 5 -^2 — — 7T, x — j -^3 



d(y,z) 



d(z, x) 



d(x,y) 







L L + — e Vect(X) 



corresponding vector field on the extended phase space X. 

Theorem 6 Vector field L e Vect(X) is a line field of the 3-form doj^ 2 \ i.e. 

i~ L (dJ 2) ) = 0. 
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Proof. Since 

duj^ = dx A dy A dz — dH\ A dH 2 A dt 

= dx A dy A dz — L 3 dx Ady A dt + L 2 dx A dz A dt — L x dy Adz A dt, 

direct calculation (cf. ||) shows that 

il{duj^) = duj^ 2 \L, ., . ) = L\dy A dz — L 2 dx A dz + L 3 dx A dy 

— L 3 (Lidy A dt — L 2 dx A dt) + L 2 {L\dz A dt — L 3 dx A dt) 

— L\(L 2 dz A dt — L 3 dy A dt) — L 3 dx Ady + L 2 dx Adz — L\dy A dz = 0. 

□ 

As in the case of Hamiltonian mechanics, this theorem has important corollary (cf. 0). 
Namely, let c be a closed 2-chain in X, g l be Nambu-Hamilton phase flow and 3-chain 
J*c = {g T c, < r < t} be a trace of the chain c under the isotopy g T . 

Corollary 




Proof. Using the formula <9(J*c) = c — g'(c), Stokes theorem and Theorem 6, we get 

[cuM- [ cu^= [ du®=0. 

Jc -V(c) JJ f c 

□ 

Next we formulate the principle of least action. Recall that in Hamiltonian mechanics 
this principle states (see, e.g., ||) that classical trajectory — integral curve 7 of Hamilton's 
phase flow with initial and final points (po, qo,to) and (pi, q±,ti), is an extremal of the action 
functional 

A(j) — [ (pdq - Hdt) 

in the class of all paths connecting initial and final points in given n-dimensional subspaces 
(if: = t ,q = g ) and (t = t\, q = qi) in the extended phase space. 
Definition 5 The functional 

A(C 2 ) = f 

—integral of the action form over 2-chains in the extended phase space X, is called action 
functional for Nambu mechanics. 
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Let 7 be a closed 1-chain in X. Define a 2-chain r in X as a trace of 7 under the isotopy 
9 f , to < t < ti, so that dT = j to —7^, where 7 to = 7, j tl = g' 1- ' ^). The following theorem 
states the principle of least action for Nambu mechanics. 

Theorem 7 A 2-chain T is an extremal of the action A(C) in the class of all 2-chains 
C whose boundaries — 1-chains c to and c tl — have the same projections onto the yz-planes as 
a given 1-chains 7t and 7^. 

Proof. As in the case of Hamiltonian mechanics, this statement follows from the previous 
theorem (cf. 0). 

Here is another proof, based on direct calculation. Assume that T admits a parametriza- 
tion x = x(s, t), y = y(s, t), z = z(s, t), where < s < 1, t < t < t\, such that it time-slices 
are closed curves, i.e. for all t functions x,y,z are periodic in s with period 1. Variations 
5x, 6y, 5z satisfy conditions Sy(s,0) = 5y(s, 1) = 5z(s, 0) = 5z(s, 1) = for all < s < 1. 
We have explicitly 

J Jn ds dt dt ds dx ds dy ds dz ds 

where II = {0 < s < 1, t < t < ti} . Using Stokes theorem and properties of 5y and 5z, we 
obtain the following expression for the variation of the action 

XA(r\ f f I ( i dz d(H u H 2 ) dy dy dijhjty dz 

5A{T) =I L { ( %-^ir%- { ¥-^^r% )5x 

dx d(H l ,H 2 ) ^dz dz , - 

1 { dt d(y,z) *ds { dt d(x,y) 1 ds } V 



d(z, 


x) 


d(H u 




d(x, 


y) 


d(H u 




9(y, 


*) 



dy d(H u H 2 ) dx dx 

+ ( %-^^r%-%-^ir } ^ } ^' 

which shows that SA(T) = for all admissible variations 5x, 5y, 5z, if T consists of integral 
curves of Nambu-Hamilton phase flow. □ 

Remark 8 The converse statement is also true: extremals of action functional are "world- 
sheets" x(s,t),y(s,t),z(s,t) consisting of families of integrals curves of Nambu-Hamilton 
phase flow parametrized by s. Indeed, these extrema are characterized by condition the the 
cross product of the following vectors in 1R 3 

(§ - {Hh H 2 , x}, I - { H h H 2 , y } , § - [Hh H 2 , 4) and (|, |, |) 

is zero, which means that there exists a function a(s,t) such that 

dx dx d(H 1 ,H 2 ) 

_ ———— 



^L-a(s t)^- = ^lllll 

dt ' ds d(z,x) 

dz_ , ,dz d(H u H 2 ) 
dt a[S ' } ds d(x,y) • 

Changing the parametrization s i— » s' = s + a(s,t)t, t \— > if = t, we see that this system 
reduces to the family of Nambu-Hamilton equations parametrized by s. 

Remark 9 Comparison between principles of least action for Nambu and Hamiltonian 
mechanics shows that "configuration space" of Nambu mechanics of order 3 constitute "two- 
third" of a phase space. 

It is straightforward to generalize presented results for the case of Nambu bracket of order 
n. The analog of Poincare -Cartan integral invariant is defined as the following n — 1-form 

u/™- 1 ) = xidx 2 A ... A dx n - HidH 2 A ... A dH n _ x A dt, 

and the action functional is given by 

JCn—l 

and is defined on the n — 1-chains in the extended phase space. In its formulation admis- 
sible variations are those which do not change projections of the boundary <9C n _i on the 
x 2 x 3 . . . x n -hyperplanes; in this case the "share" of "configuration space" in a phase space is 
1 - 1/n. 

Remark 10 This construction of action form and action functional is somewhat similar 
to the construction of cyclic cocycles in Connes' approach to the non- commutative differential 
geometry 0. 

4 Quantization 

There exist several different (and, in a certain sense, equivalent) points of view on quantiza- 
tion problem. 

One is based on the approach which uses a deformation theory of associative algebras. It 
considers quantization as a deformation of a (commutative) algebra of classical observables 
on the phase space in the "direction" defined by a given Poisson (or symplectic) structure 

0- 
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Namely, let X be a Poisson manifold with the Poisson bracket { , } and algebra of classical 
observables A = C°°(X). One-parameter family {Ah} of associative algebras is called a 
quantization of a commutative algebra A of classical observables, if following conditions are 
satisfied. 

1. Algebra A is included into this family 

Ah\h=o = A, 

(or, in a formal algebraic category, A = Ah/hAh). 

2. All algebras Ah are isomorphic to A (or rather to ^4 [[/?<]] in formal algebraic category) 
as linear spaces. 

3. Denoting by *h (associative) product in A h (^-product in the terminology of |]I2"1), one 
has the expansion 

/i*h/ 2 = /i/ 2 + |c 2 (/ 1 ,/ 2 ) + o(/ i 2 ), 

with the property 

C 2 (/ 1) / 2 )-C 2 (/ 2) / 1 ) = {/ 1 ,/ 2 }. 
The last property is often referred as a correspondence principle 

{fx, f'2} = lim \(f x * h f 2 - f 2 *h fi) 

between classical and quantum mechanics. 

Deformation theory approach can be carried out explicitly in the case X = R 2 (or R 2n ) 
with canonical Poisson bracket and yields celebrated Hermann Weyl's quantization scheme. 
Corresponding ^-product *h — a map from A £g> A into A, is given by a composition of a 
usual commutative point-wise product on A with the following bilinear pseudo-differential 
operator 

exp(^{ , }) = exp(^{^- A ^-}) :A®A^A®A. 

One might use the same approach towards the quantization of Nambu mechanics. Namely, 
suppose that we try to include "usual" ternary product — a map from A <g> A <g> A into A, 
given by a point- wise multiplication of classical observables, into the "new" ternary operation 
( , , )h which depends on parameter h and satisfies certain natural properties (analogous to 
associativity condition in the binary case) and a correspondence principle 

{/i,/a,/s} = Bm^Alt(/ 1 ,/ 2 ,/ 3 ) fcl 
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where Alt denotes complete anti-symmetrization with respect to the arguments /i, /2, f%. 

We can easily generalize Weyl's formula by letting the new ternary operation (,,,)& 
to be the composition of a point-wise ternary product and the following trilinear pseudo- 
differential operator 

exp(^{, , }) = e X p{^{^A^A^}):A9A9A^A»A9A. 

This formula generalizes verbatim for canonical Nambu bracket of order n. 
Remark 11 One should note the following. 

(i) Contrary to the usual binary case, point-wise ternary product of classical observables 
looks less natural. 

(ii) "Natural" constraints for ternary operations generalizing usual associativity con- 
straint are, apparently, not well understood. 

Nevertheless, proposed deformation has certain "appeal" and should be analyzed further. 
It might lead to the (partial) answer to a problem in (ii). 

Remark 12 Recently R. Lawrence |H| proposed a system of axioms for linear spaces 
with n-aiy operations (n-algebras using her terminology). It looks that "deformed" product 
( , , )h, as well as a ternary point-wise product, do not satisfy them. However, one should 
have in mind that approach in |H| generalizes certain combinatorial properties of algebraic 
operations (notably the Stasheff polyhedron), whereas our approach is based on "dynamics". 

Another approach uses Feynmann's path integral formulation of quantum mechanics. For 
a quantum particle, described classically by a Hamiltonian system, it gives the probability 
amplitude of the transition from a state \qo) at time t = to into state \qi) at time t = t% as 
a functional integral of the exponential of the classical action 

exp{-A(7)} 

over all "histories" 7 with respect to the "Liouville measure" 

i 

—DpDq 
h 

in the functional space of all "histories". 

Principle of the least action for Nambu mechanics can be used to formulate just similar 
rule. However, special form of action principle (see Theorem 7) requires that quantum 
states \y(s),z(s)) should be parametrized by loops y(s),z(s) rather then by points in a 
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"configuration" space. This natural appearance of loops looks quite appealing. It suggests 
that one might still have particles as a point-like objects in classical picture and dynamic of 
loops (or, generally, n — 2-closed chains) in quantum picture. We are planning to investigate 
this intriguing opportunity elsewhere. 

Finally, we present yet another approach to quantization a canonical formalism. It is 
based on Heisenberg commutation relations, which for the phase space X = H 2 with canon- 
ical Poisson bracket look the following (in a complex form) 

[a, a'] = aa) — a) a = I, 

where operators a', a act in a linear space of quantum states. They have the following 
realization in a space H2 with the basis {|n)} n >o parametrized by non-negative rational 
integers: 

a\n) = \/n|n — 1), a)\n) = \/n + l\n + 1). 

The vector |0) plays the role of a vacuum state and operators a', a are called creation- 
annihilation operators. 

Being one of the fundamental principles of quantum mechanics, Heisenberg commuta- 
tion relations have remarkable mathematical properties. In particular, one has celebrated 
Stone-von Newmann theorem that all irreducible representations of Heisenberg commuta- 
tion relations are unitary equivalent. Application of Stone-von Newmann theorem to the 
linear canonical transformations gives a certain projective representation of the symplectic 
group, which was shown by Andre Weil [nj] to play a fundamental role in the arithmetics of 
quadratic forms and quadratic reciprocity for number fields. 

In (]]] Nambu proposed the following generalization of Heisenberg commutation relation 

[A 1: A 2 , A 3 ] = A X A 2 A Z - A X A^A 2 + A Z A X A 2 - A 3 A 2 Ai + A 2 A 3 A 1 - A 2 A X A Z = cl, (8) 

where Ai,A 2 , A3 are linear operators, J is a unit operator and c is a constant. We will call 
(§) Nambu- Heisenberg relation. 

Remark 13 Although Nambu-Heisenberg relation satisfies correspondence principle 

£[,,]->{,,}, 

a priori it is not clear that observables in quantum Nambu mechanics should be realized by 
linear operators. Therefore the form ([8]) might look rather ad hoc. It was already pointed 
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by Nambu |jj that, contrary to the usual commutator [ , ], "ternary commutator" [ , , ] 
is not a derivation with respect to the operator product, which poses certain problems in 
formulating quantum dynamics (see M for detailed discussion). From our point of view, 
it looks like there are no reasons to impose such derivation property with respect to the 
operator product. "Triple commutator" [ , , ], should be considered to define a Nambu- 
Lie "gebra" structure on quantum observables. Its property of being a derivation of this 
structure (see Definition 3) will imply consistency of quantum Nambu dynamics. 

Instead of continued this discussion (which we are planning to do elsewhere), here we 
will only present explicit realization of Nambu-Heisenberg relation @. 

Let Q[p] be a quadratic number field with 1 + p + p 2 = 0, i.e. 

-1 + 7=3 
P= 2 ' 

and let Z[p] be a ring of algebraic integers in Q[p], i.e. 

uj = mi + m 2P G Z[p], mi, m 2 G Z. 

Denote by H 3 linear space with the basis {\u))} parametrized by Z[p]. Direct calculation 
proves the following result. 

Theorem 8 Nambu-Heisenberg relation with c = p — p = ^—3 admits the following 
representation 

Ai\u) = (u + 1 + p)\u + 1), A 2 \uj) = (u + p)\u + p), A 3 \u) = uj\u + p 2 ) 
in the space H 3 . 

The following arguments explain why such type of reprentations are possible. 
Using the analogy with n = 2 case, consider representation (|j) of the form 

A 1 \u) = f 1 (u)\u + 1), A 2 \u) = f 2 {u)\u + p), A 3 \u) = f 3 (u)\u + p 2 ). 

Because of 1 + p + p 2 = 0, vectors \u) are eigen-vectors for all possible triple products 
Ai 1 Ai 2 Ai 3 of Ai,A 2 ,A 3 and Nambu-Heisenberg relation (||) reduces to functional equation 
for the uknowns fi(u>), f 2(^)7 fsi 1 ^)- Theorem 8 exhibits one of its solutions. We can obtain 
other solutions using the following trick. Equation @ is invariant under the similarity 
transformation 

A i »A i = U- 1 A i U, * = 1,2,3, 
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where U : H 3 i— > H 3 is an invertible linear operator. Choose U to be diagonal in the standard 
basis of H 3 , i.e. 

U\u) = u(uj)\uj), u(u) 7^ 0, u e Z[p]. 
Then operators A { are represented by the same formulas as A^s with 

j\(uj) = u' 1 ^ + f 2 (u>) = u~ l (u) + p)u(uj) f 2 (uj) , f 3 (u;) = u~ 1 (uj + p)u(u>) f 3 (u;) . 

Now assume that A\ and A 2 commute and /i,/2 are non-zero for all u G Z[p]. Then one 
can choose u{uj) such that fiiyj) = f2(u) = 1- Indeed, u{uj) can be found as a solution of 
the following compatible system (because of the commutativity of A± and A 2 ) of difference 
equations 

u(u + l) = fi(u)u(u), 
u(u + p) = f 2 (u)u(u). 

(Provided certain conditions at "infinity" for the functions f\ and f 2 are satisfied, u{uj) 
can be given as semi-infinite product over the lattice Z[p].). Therefore, original functional 
equation for /j's reduces to the following simple equation 

/ 3 (W + 1)-/3H = C, 

which can be easily solved. 

"Gauging back" these solutions with /i = f 2 — 1 we obtain other solutions. 

Solution presented in Theorem 8 does not belong to this class: operators A 1 and A 2 (as 
well as other pairs) do not commute and, since vectors |p 2 ), |1 + p 2 ) and |0) annihilate, 
operators A ± ,A 2 and A 3 , corresponding functions /j have zeros. This implies that in this 
case corresponding invertible operator U does not exist. This shows the special role played 
by representation in Theorem 8: although there is no vacuum vector in a conventional sence, 
there exists triple of vectors 

|p 2 >, + A |0> 

with "similar" properties. 

Remark 14 It should be noted that this realization of Nambu-Heisenberg relation must 
be supplemented with the analog of Hermitian anti-involution, which is this case should be 
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an operation of order 3 (since we a dealing with a triple instead of a pair). We are planning 
to address this issue elsewhere. 

Finally, consider briefly the general case of Nambu mechanics of order n. One can 
postulate the following form of Nambu- Heisenberg relation 

J2 (-iy^A a{1) ...A a{n) = cI. (9) 

aaSymrn(n) 

Assuming that n is a prime integer greater then 2, consider representation of (Q) in a linear 
space H n with the basis {\uj)} parametrized by algebraic integers Z[p] in the cyclotomic field 
of n-th root of unity, i.e. p n = 1, or 

1 + P + ... + P n - 1 = 0. 

Assuming that operators Aj are represented in a similar way: 

Ai\u) = fi(uj)\uj + p l ~ l ) , i = l,...,n, 

one might try to repeat the same trick: "gauge out" the coefficients fx, ... , f n -\. Namely, 
assume that f^s are such that the following system of equations 

u{l) + p 1 ^ 1 ) = fi(u>)u(u>), i = 1, ... n — 1 

is compatible and has a solution which gives rise to an invertible operator U. In this case 
commuting operators Aj = U~ 1 A i U for i — 1, . . . , n— 1, will represent n— 1 basic translations 
(generators) of the lattice Z[p] and, therefore. However, analyzing relation (|9]) for general 
n we see that all its terms contain at least one pair of A^s with i < n — 1 as nearest 
neighbors. These terms will appear twice in (^) and commutativity of Aj's implies that they 
will mutually cancel each other! 

This shows that in the higher order case Nambu-Heisenberg relation does not admit 
simple solutions with n— 1 commuting operators A^s. However, we suspect (and are planning 
to investigate it) that it has solution similar to that in Theorem 8. 

5 Conclusion 

Though this paper poses more questions that provides the answers, we feel that the subject of 
higher order algebraic operations might be relevant for future development of mathematical 
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structures related to physical problems. We suspect that it might clarify certain problems 
related to the generalizations of the integrability concept (Yang-Baxter equation, Poisson- 
Lie groups, quantum groups) for the higher dimensional case (Zamolodchikov's tetrahedron 
equations, Frenkel-Moore solutions, 2-categories) . Ending on highly speculative note, one 
might suggest that perhaps these higher order structures can be also relevant in the arith- 
metics of forms of higher degrees, higher reciprocity laws in number theory, etc. 
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